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The effects of variations in the components of a beam waveguide are con- 
sidered. These variations statistically cause the Gaussian beam spot size 
of the light propagating down the waveguide to grow and cause the beam 
center to oscillate about the waveguide axis with ever-increasing amplitude. 
Random variations in lens focal length and spacing and random lateral 
lens displacements are considered. It is also shown how random variations 
in focal length and spacing can be included in the published analyses for 
short random bends in the waveguide axis. 

When the number of lenses is large, it is shown that the beam displace- 
ment and beam spot size both grow exponentially with distance. 

As an example, a confocal waveguide with lenses spaced one meter apart 
and built to somewhat optimistic tolerances will require a beam redirector 
every 2.5 kilometers to prevent the beam oscillations from exceeding an rms 
value of 2 millimeters. 



I. INTRODUCTION 

A long sequence of spaced lenses is of considerable interest for optical 
connnunications. It is known that the diffraction losses in such an 
optical beam waveguide can be kept very small for moderate size 
lenses. 12 This means that if a transmission line is made of identical 
low-loss lenses, spaced identically along a straight line with each lens 
centered on this straight line, there is a mode of propagation which is 
low loss. However, if there are imperfections in the transmission line, 
the light beam will begin to wander from the axis or the beam size will 
grow and the beam will eventually strike the edge of the lens and be 
lost. Since the diffraction loss of the beam in a perfect line can be kept 
very small, it is the line imperfections, the line axis curvature, and the 
scattering and absorption at each lens which will primarily determine the 

451 



452 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1966 

optical loss. Gas lenses have been considered for reducing the scattering 
and absorption losses. 3 

Rowe 4 and Hirano and Fukatsu 5 have shown how the beam position 
is affected by random lateral lens displacements. Berreman, 6 Marcuse, 7 
and Unger 8 have considered correlated lateral lens displacements in the 
form of bends. All of these analyses have assumed perfect lenses and 
perfect spacing and have shown the growth of the beam displacement 
due to lateral lens displacements only. It is the purpose of this paper to 
show how the previously obtained results are altered when the lens 
focal lengths and lens spacings have random variations. 

In this paper, we shall consider the statistical effect of random varia- 
tions in lens focal length and lens spacing and random lateral lens dis- 
placements. We shall also consider random bends whose correlation 
length is much smaller than the total line length. It is shown how the 
various line imperfections couple to one another and cause the beam 
deviation from the axis to grow. The growth of the spot size of a Gaus- 
sian beam is also considered. 

It is shown here that the random variations in / and L cause the rms 
expected value of beam displacement and the rms expected value of the 
beam spot size to grow exponentially with distance when the number of 
lenses is large. For bends, the variations in / and L cause an exponential 
increase in the average allowed bending radius of the guide when the 
number of lenses is large. In contrast, when / and L are perfect these 
effects grow more slowly with distance, and increase only as the square 
root of the number of lenses. 

We shall use geometric optics since it is known that in the paraxial 
approximation the center of a Gaussian beam in a beam waveguide 
behaves as a ray. 9 The geometric optics analysis is extended to find the 
behavior of the beam spot size by replacing the Gaussian beam by its 
equivalent ray packet. 10 

II. GENERAL PROBLEM FORMULATION 

We shall consider the problem in two dimensions only for simplicity. 
It has been shown that the three-dimensional problem can be split 
into two independent two-dimensional problems. 5 For aberration-free 
lenses, the motion of the beam in one transverse dimension is dependent 
only on the initial conditions and lens displacements in the same trans- 
verse direction. 

Consider the transmission line shown in Fig. 1. We define r„ and r„ 
as the position and slope of the ray just to the right of the nth lens. The 
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Fig. 1 — Beam waveguide notation. 

ray position and slope are measured with respect to a straight line 
which is the nominal transmission line axis. The spacing between lenses 
is labeled L n and the convergence of the lenses as C n where 

Cn = 1/fn 

and /„ is the focal length of the nth lens. The lateral distance between 
the center of the nth lens and the reference line is s„ . The displacement 
s„ is positive if the lens center is above the reference line. 
Using this notation we can write 

r„ = r„_i + LnTn-i (1) 

r„' = -CV,,-! + (l-C n L„)r„_/ + C n s n . (2) 

If we define 



Rn = 



M n = 



and 



V H = 



1 L n 



C n S n 



we can write (1) and (2) in matrix form as 

R n = MJtn-1 + Vn ■ 



(3) 



This relates the ray position and slope at the n plane to the ray position 
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and slope at the n — 1 plane. We shall be interested in the rms expected 
value of the output beam displacement and hence in the square of the 
output beam slope and position. We shall, therefore, square the matrix 
(3). To do this, we take the Kronecker product 11 of each side of (3) with 
itself. 

En X Rn = (M n X M„)(R n -l X Rn-l) + V n X V n 

v (4) 

+ (MnR,,-!) X V n + V n X {M n R n ^) . 

Now take the expected value of (4) 

(Rn X Rn) = ((Mn X M„)(R n -i X Rn-*)) + (Vn X V n ) 

+ ((MnRn-l) X Vn) + (V n X (M n R n ^)). 

The expected value of a matrix is the matrix of the expected values. 

We will now state our statistical assumptions. We assume small 
variations about L and C 

L n = L(l +Z„) 

C n = C(l +c„). 

(In) = (Cn) = (Sn) = 
(Uk) = (CnCk) = n 7* Ic 

(InCk) = (InSk) = (CnSk) = all 71 

I 2\ 2 

(C„ ) = Oc 

(L 2 ) = <r L 2 

(s n 2 ) = y\ 

Hence, l n , c„ , and s„ are mutually independent random variables of 
zero mean. There is no correlation between the variations in spacing, 
the variations in focal length, and the variations in lateral displacement. 
The variations in L and C are completely random with no correlation 
between adjacent L's and adjacent C's. For the first sections of the 
paper we shall consider the lateral lens displacements to also be com- 
pletely random with no correlation between adjacent displacements. 
In a later section it will be shown how correlation of the lens displace- 
ments in the form of random bends in the waveguide axis can be included. 
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III. A TRANSMISSION LINE WITH RANDOM LATERAL LENS DISPLACEMENTS 

In this section we shall consider only random lateral lens displace- 
ments. Hence, we impose the additional statistical restriction that 

(s k s n ) = n ?* k. 

The lateral displacement of any lens is unaffected by the lateral dis- 
placements of any other lens. 

If we repeatedly substitute (3) into the last two matrices of (5), 
we see that they contain elements of the general form 

(G(L k , Lfc+i , • • • L» , Ck , Cfc+i , • • • C„)s*_is n ), k ^ n 

and 

<F(Li , U , • • • L„ , d , C 2 , • • • C„)R s n ) 

where G( ) and F( ) are some functions. In view of the statistical 
assumptions, these can be written as 

(G(L k , L k+1 ■■■ L n , C k , C k+1 ,■■■ C n )>(s fc _ lS „), k ^ n 

and 

(F(ln , U , • • • L n , Cx , C t|1 . ■ • C n ))R (s n ) 

which are zero since (s fc _is„) = for k ^ n and (s„) = 0. Hence, the 
last two matrices of (5) are both zero. 

When there is some correlation between lateral lens displacements, 
i.e., a wavy transmission line axis, these two matrices are not zero. It 
is through these matrices that the correlation will enter. 

Also, by using the above statistical assumptions M„ is independent of 
#„_! and (M„ X M n ) and (V n X V n ) are not functions of n. Equation 
(5) can, therefore, be written as 

(R n X R n ) = (M n X M n ) n R X Ro + E (M n X M n f (V n X V n ) (6) 

where 7? is the matrix of the initial ray slope and position. 

The Kronecker products in (6) are either 4 X 4 or 4 X 1 matrices. 
These can be reduced to 3 X 3 and 3X1 matrices by combining the 
two redundant terms." For clarity and ease of computation we shall 
reduce the matrices and write them out explicitly below. We have 
assumed here that <r L 2 « 1 and a c 2 « 1 hence, we have neglected 
<T L 2 a c 2 compared to a L 2 or cr c 2 . 



456 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1966 

(R n X Rn) = (r n r n ') 
Ar?) . 



R X R = 



To' 



r r 



r 



(V n X V n ) = 







(M„ X M n ) = 



1 

-c 



|_C 2 (1 + cr c 2 )2/ 2 J 
2L L 2 (l + err, 2 ) 

1 - 2LC L - L 2 C(1 + (Tz, 2 ) 



C 2 (1 + (T C 2 ) -2C + 2LC 2 1-2LC + L 2 C 2 

(1 + <rc 2 ) (1 + *l +<rc 2 ) 



3.1 27ie Characteristic Roots of (M n X M„> 

To evaluate (6) will require the raising of (M n X M n ) to the kth 
power. To do this it will be necessary to find the characteristic roots of 
(M n X M n ). The characteristic roots of (M n X M„) can be found from 
the equation 

| (Mn X M n ) - IX | = 

where X is the characteristic root and J is the unity matrix. This leads to 
the following cubic equation for X : 



(7) 



X 3 - X 2 [3 - 4LC + L 2 C 2 + L 2 C 2 (cr L 2 + a c 2 )] 

- X [-3 + 4LC - L 2 C 2 + L 2 C 2 (a L 2 + <r c 2 )] -1=0. 

We assume in (7) that cr L 2 and ac are very small quantities, hence 
terms of higher power than 2 in a L and a c are neglected. Since a L 2 
and a c are assumed very small it is reasonable to assume that the roots 
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of (7) are very near the roots for the perfect transmission when a L = 
a c = 0. For the perfect transmission line, the roots are 

X = 1, e , e , 

where = cos -1 (1 — LC/2). We therefore write the three roots of (7) 

as 

Xi = 1 + ft 

X 2 = e 2,9 (l + 32) 

X 3 = <T 2,9 (1 + g 3 ) 
where 

l(?i|, l&li I 03 1 « 1. 
For the case of LC ^ 2, i.e., a nonconfocal system, (7) gives 

2LC / 2 , 2\ 
0i = a _ Tf t to"£ + <rc ) 

— LC , 2 , 2v 



If we define 



then 



92 - 93 = ^ _ LG (<rz. + <rc ). 



2jL(7 / 2 1 2\ 





Xi = 


1 + a 






X 2 = 


e UB (l - 


■i> 




X 3 = 


e-™ (l 


-!)• 


For the confocal case 


LC = 


= 2, the 1 


oots of (7) are 




Xi = 


1 + a 






x 2 = 


-1 + a 






X3 = 


-1. 





These two sets of solutions of (7) are valid so long as a « 1. This 
means a c « 1, ol « 1, and LC is not near 4. 



(Xi 


- x 2 


)(Xi- X.) 




+ 


U- 


-Xi/)U - 


-X 3 I) 


(X. 


- Xi)(X 2 - 


- X,) 


_l_ 


(A- 


-Xx/)U - 


-X 2 7) 
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3.2 Sylvesters Theorem 

For raising the matrix (M n X M n ) to the /eth power, it is helpful to 
use Sylvesters Theorem. 12 If Xi , X 2 , and X 3 are the characteristic roots 
of the matrix A then 

A k = U - MI) (A - \ 3 I) ,k 



A2 



X* 

1 (X s - Xi)(X 8 - Xi) 3 

where I is the unity matrix. 

In (6) we shall be interested in the case where n is a very large num- 
ber, i.e., many lenses in the transmission line. The difference between 
the matrix (M„ X M n ) and M X M where M is the perfect transmission 
line matrix is very small (terms of the order of a c 2 and a L 2 ). For the 
nonconfocal case, it will therefore only be necessary to consider the 
deviations from a perfect transmission line in Xi , Xa , and X 3 . In the 
coefficients of Xi*, X 2 fc , and X 3 * we can assume <r L 2 = a c = 0. 

Hence, for the nonconfocal case, we can write 

(M X M f = (M X M ~ e "° I)(M XM ~ e " Ul) r" 

\m. n a nin) (1 — e 2ie )(l — e~ 2ie ) 

(M X M - I)(M X M - e**I) **>m (h) 

+ (g2ie _ 1 )( e 2.-9 _ e -2f9) <> 

(MXM - I)(M XM - e 2i0 I) -tm-w 
+ ( e *«- l)(e-™ - e™) 

Here we have used 1, e 2 ' 9 , and e~ 2,e for the roots of M X M and 

Xi* = (1 + a) k « e ka 

X 2 k = e 2 '™ (l - ?Y « e 2, 'V fc " ' 2 



\ V 



X 3 fc = <T 2rt9 (l -?] ««"**'«" ta/1 



In (6) we shall be interested in only the first element of (R n X i2») 
which is (r 2 ). To calculate this we will only need the first row of ele- 
ments of (M n X M n f. These can be found from (8). The 1,1 element for 
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the nonconfocal case is 

ka f4 sin" 6 - 2LC cos 6~] 
|_ 4 sin- 6 J 

, tikt -kail [ -2LC cos 6 1 
e e L(e 2 '' 9 - 1)2.; sin 20 J 

-■like -kaii [~ -2LCcosfl 1 

"'" 6 LCe -21 * - 1) ( -2j sin 20) J ' 

After some simplification 

611 = T=~lc [e * fl + e ~* a ' 2 cos (2A; ~ 1)0] - 

Similarly, we can write b n and 613 for the nonconfocal case as 

'» = r^rc [/" + e ~"°" (Vt!F sin 2W - cos M )] 

For the confocal case, LC = 2, we are close to a degenerate case 
where two characteristic roots are close to being equal. If we retain 
terms to no higher power than 2 in <r L 2 and oc in (M n X M n f we can 
write 

{Mn x ^ m (MXM + IUMXM + I) e , 

({M n X M n ) - (1 + a)/)(<M» X M n ) + /) (__! )»,-*. (g) 
2a 

((J/, X il/„) - (1 + a)I)i(M H X M») + (1 - a)/) ,_ )k 
+ 2a * 

We have retained the a* and a c 2 in the coefficients of the last two terms 
since vl and <r c 2 appear in the denominators. 

We again are interested only in the first row of (M„ X M„)'. For the 
confocal case, the elements of the first row are 

j. ka 

i r I ka 2( — 1)/ 2 -ka , 2 \ 

Dm = L\ e ((fee + a L ) 

, L I to 2(— 1) / 2-Jfca . 2\ 

3 = ~2 
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3.3 Calculation of (r„ 2 ) 

For the nonconfocal case using (6) and the values for the elements 
in the first row of (M„ X M n ) k , we can calculate (r n a ). 

(r« a ) - :r ^ w L 2 [e na + e - (no/2) cos (2n - 1)0] 

+ r Lr ' (V + <T (na/2) (a/ 4 ~£ C s i n 2n0 - cos 2nfl)l 
+ &£. y _ ,-*■» C o S 2n0] 



+ W Z (e ta - e~ (ia/2) cos 2*0)1. 
t=o J 



In the 3,1 position of (V n X F n ) we have neglected the a a as compared 
to 1. 

The summations can be evaluated as 

n— 1 na i na i 

V-> ka __ c - l^e — 1 
6o 6° — 1 ~ a 



n-l 

E 

fc=0 



E e" (fca/2) cos 2/c0 

_ 1 - <T (a/2) cos 20 - <r ( " a/2) cos 2n0 + e - 1(n+1)a/21 cos 2(n - 1)0 
1 + e-° - 2e-° /2 cos 20 

1 . e- (nal2) sin (2n - 1)0 
W - H . 

2 T 2 sin 

We have used the facts that a <JC 1 and n»l. 

The expected value of the square of the output beam position for the 
nonconfocal case is therefore, 

<r„ 2 > = r= ^ jroV + e- (na!2) cos (2» - 1)0] 

+ roLro' [V + «~ (M/1) (a /* ~^° sin 2n0 - cos 2ndX\ 



T 2 . 2 

Xj 7"o r no — (na/2) nl 

H — p^ Le — e cos 2n0] 

-L/C 

, rn 2 IV - 1 1 e~ {nal2) sin (2n - 1)0] 
+ LC2/ |_ _ a 2 2-slnT - J 



(10) 
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where 

2LC / 2 - 2\ 
a = T — w (a L + a c ). 

For confocal spacing, LC = 2, the square of the expected value of 
r n is 

A. 2 \ ». 2 o no I „. r-' „ na 2(— 1) / 2 -na , 2\ 

\r n ) = r e + r Lr e 1 (o- c e + <r L ) 

I L To „ a 2( — 1) , 2 -no . 2\ / 1 i \ 

H 2~ a — c +^) C 11 ) 

[„no i i 2 —na / ■§ \n 2~| 

,_-! _ 1 aoe_ + (dO«t] 
a 2 a a J 

where, for LC = 2, a = 2(«r t 2 + (r c 2 ). 

It is of interest to consider how close to LC = 2 one must be to have 
(11) hold rather than (10). If we retain terms to only the first power 
in a, it can be shown that (11) is valid when 

\LC - 2 | < a L 2 + <rc 2 - 

If 

| LC - 2 | > * L 2 + ac 2 

then (10) holds. Since a L 2 and a 2 will be of the order of 10 -4 we must 
be very close to LC = 2 for (11) to hold. 

If the lenses and spacing are perfect so that a = 0, the first three terms 
of (10) and of (11) give the square of the output beam position due to 
the input beam slope and position. The last term gives the increased 
displacement due to random lateral lens displacements. Both parts 
agree with Hirano and Fukatsu 5 when a = 0. 

The random errors in focal length and spacing cause an exponential 
increase in the expected value of the square of the output beam dis- 
placement. This can be seen more clearly for the case where n is very 
large. If na > 2, then (10) and (11) reduce to the same result. In this 
case, 

< r " ) ^4^LcL r ° +roLr ° + TcJ + 4^lc4~^-J- (12) 

IV. TRANSMISSION LINE WITH RANDOM BENDS 

We will now assume the axis of the beam waveguide is bent so that 
there is some correlation between adjacent lateral lens displacements. 
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As noted in Section III this correlation will appear in (5) in the last 
two terms, ((M«R n -i) X V») and (V„ X (M„fl»_i)>. 

It is shown in the Appendix that these two terms to first order do 
not contain <r L 2 or a c 2 , i.e., they are not affected by the random varia- 
tions in focal length and spacing. It is also shown that if the axis of the 
guide is composed of a series of uncorrelated bends whose average bend 
length is much smaller than the total length of the transmission line, 
these two terms are not functions of n. This type of bending might 
typically be the case for a very long transmission line laid to follow the 
gentle bends of the terrain. 

It was shown in Section III that to first order (V„ X V„) also has 
these properties, i.e., it is not a function of <t l 2 , <tc , or n. Because of this 
similarity, the matrix ((MX-,) X V n ) + <F„ X (M n Rn-i)) can be 
considered as an added part of (V„ X V„) and can be carried through 
the analysis in this manner. Hence, the random errors in focal length 
and spacing affect the beam displacement due to short uncorrelated 
bends in the same way they affected the beam displacement due to 
random lens displacements. 

From (10) or (11) for "a" small we can show how ol and ac couple 
to the random displacements by writing 

,2, / «% e — 1 

(r n )= (r n ) a=0 —^-. 

In this expression, <r„ 2 ) =o is the expected value of the square of the 
beam displacement due to random lens displacements when a = 0. 
Because of the similarity pointed out above, the beam displacement 
due to short random bends is also multiplied by (e" a — \)/na to account 
for the focal length and spacing errors. Let us assume a transmission 
line axis is specified which fits the conditions, i.e., it is composed of a 
series of uncorrelated bends whose bend length is much shorter than the 
total line length. From this we can calculate (r„ 2 ) assuming L and C are 
perfect. This has been done for some specific cases by Marcuse, Berre- 
man, 6 and Unger. 8 To include random imperfections of L and C if a is 
small we multiply this result by 

na 1 

e — 1 



na 



This analysis does not hold if the correlation extends along the entire 
line (for example a serpentine bend) or if the correlation extends over a 
large portion of the line. 
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V. STATISTICAL GROWTH OF BEAM SPOT SIZE 

We have been concerned thus far with the behavior of light rays in an 
imperfect transmission line. Our primary concern, however, is the 
behavior of Gaussian light beams rather than light rays. It has previously 
been shown that in the paraxial approximation the center of a Gaussian 
light beam does behave like a ray. 9 We can regard r n and r„' therefore, 
as the position and slope of the center of a Gaussian beam at lens n and 
r and n' as the initial conditions of the beam center. We lack infor- 
mation on the effects of the transmission line imperfections on the size 
of the beam. 

Using the complex beam parameter law of Kogelnik 13 it would be 
possible to find the statistical growth of the spot size due to the line 
imperfections. It will be simpler, however, to use Steier's ray packet 
equivalent to the Gaussian beam. 10 This approach conveniently uses 
the already derived statistical behavior of the light rays to find the beam 
size behavior. 

Just to the right of a lens, the ray packet equivalent of the normal 
Gaussian mode of the transmission line is 



i L ■ 
?*o = w cos ip + r. — sin (p 

k,wo 



-2 . 

ro = i — sin <p 
kwo 



(p has all values from to 2w (13) 



where 



wo = spot size at the beam waist = — - — r— — - 

L kC J 

lc ~T 

X = wavelength. 

If the path through the transmission line of each ray of the packet is 
found then the behavior of the equivalent Gaussian mode through the 
transmission line can be found. At any point in the transmission line, 
the envelope or the distance between the extreme rays of the ray packet 
is equal to twice the beam spot size and the curves which are perpendicu- 
lar to the average ray slope are the beam phase fronts. 

To find the effect of transmission line imperfections on the beam 
spot size, let us launch the ray packet given by (13) into the trans- 
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mission line. If we substitute the value for r and r ' from (13) into (10) 
for the nonconfocal case we find 

<>••'> = t^tc e " + A «- < "" 2 ' [cos 2 * - (2n - l)e d4) 

— sin 2<p sin (2n — 1)0] 

where tp ranges from to 2tt. We have not included the last term of (10) 
since the lateral lens displacements have no effect on the growth of the 
spot size. 

The expected value of the square of the spot size at the nth lens, 
(w n 2 ), is given by the envelope of these rays. Taking the maximum value 
of (14) as <p goes from to 2?r. 

. 2 \ 2Wo /no | -(no/2)\ 

(u>« ) = i - LC (e + e ). 

The normal mode spot size at a lens, w, is given by 

2 4w<> 2 

w T^W 

Hence, for the nonconfocal case 



where 



/ 2\ no , — (na/2) 

(«Q e + e / 15 n 

w* 2 



2LC , 2 , 2 - 



AU\j / 2 ■ 2\ 

4 — LC 



For the confocal case, LC = 2, we substitute from (13) into (11). 
Taking the maximum value as <p ranges from to 2ir. 

, 2\ -. ( 2 —na , 2"\ 

<SE*1 „IL"- + f2Lf ± ^ (16) 

If na is small so that 

e na tt 1 + na, 
then for the nonconfocal case (15) becomes 

(Wn) , , LCn , 2 - 2\ / 1<7 \ 

V =1 + 2(4-LC) ( ^ +irc) ' (17) 
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and for the confocal case (16) becomes 

H 1 = 1 + «*'. (18) 



Equations (17) and (18) agree with the results of the perturbation 
analysis of Hirano and Fukatsu. 5 

Hence, for na small, the random errors in C do not affect the spot 
size in the confocal case. However, as pointed out in Section 3.3, LC 
must be almost 2 for this to be true. If | LC — 2 | > a J + a c 2 the result 
for a nonconfocal system should be used. Since a L 2 , o- c 2 ~ 10 -4 , this is 
a very stringent requirement on LC. It is doubtful if LC can be held 
close enough to 2 to gain this advantage in reduced spot size growth. 

If na > 2 then the nonconfocal and the confocal results are very 
nearly the same and for both cases 

W) „ 1 na 
VI. SUMMARY 

The results derived here show statistically how imperfections in an 
optical transmission line affect the output beam from the transmission 
line. The imperfections cause the beam center to wander from the 
transmission line axis and cause the beam size to grow. We have con- 
sidered the errors in focal length and spacing to be random and the 
lateral lens displacements to be random or with short correlation lengths. 
For this case, statistically the size of the beam and the distance of the 
beam center from the axis grow exponentially as the number of lenses 
when there are many lenses. 

For na > 2, and random lateral lens displacements, the results can be 
summarized as follows. The beam center launched at r , r ' has an rms 
expected value of 

The beam spot size has an rms expected value of 

- H 9\ "a/2 W 



For random transmission line bends of short correlation length the 
random variations in L and C increase the expected value of the square 
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of the output beam displacement as 

. 2 X e —1/2. 

<<■„ > = -j^- (r. >„» 

where (r n 2 )a-o is the value computed assuming no variations in L and C. 

If a beam is launched into a straight line on axis with no slope its 
position is not affected by random errors in / and L and only the size 
of the beam is affected. This is obviously true since a ray through the 
center of a lens does not bend no matter what the lens focal length. 
If, however, the axis is curved or the lenses have random lateral dis- 
placements, the beam begins to wander from the axis and is now affected 
by the errors in / and L. This coupling is clearly shown in these results. 

These calculations are pertinent when n is large. This is the case of a 
transmission line with relatively closely spaced lenses which would be 
able to control the light beam around gentle bends in the terrain. 

As an example, let us consider a confocal beam waveguide with lenses 
spaced every one meter and built to the following rms expected value 
tolerances : 

(i) focal length variations — 1 per cent 
(it) spacing variations — much less than 1 per cent 
(Hi) random lateral lens displacements — 2 X 10~ 2 mm. 
These tolerances give 

a = 2 X 10" 4 

x f = 4 X lO" 10 m 2 . 

If we assume an rms output beam deviation of 2 millimeters is 
acceptable, we can go approximately 3.5 kilometers (n = 3.5 X 10 3 ) 
with the line described above. If the line is allowed to have gentle 
circular bends of an average radius of curvature of 5 kilometers and 
an average bend length 7 of 100 meters then the distance which can 
be traveled before there is an rms beam deviation of 2 millimeters is 
reduced to 2.5 kilometers (n = 2.5 X 10 3 ). This means a beam redirector 
is required every 2.5 kilometers. 

In the above example, the dominant Gaussian mode spot size at each 
lens for the perfect line is 0.45 mm rad. Because of the line imperfections, 
this grows to an rms expected value of 0.48 mm rad at n = 2.5 X 10 3 . 
This spot size growth is insignificant compared to the beam deflection. 
In general, unless LC is very small the spot size growth is not as im- 
portant as the beam deflection growth for closely spaced lenses. 

In calculating these numbers we have assumed the lenses are perfect 
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and have neglected any aberrations. Additional work is required to 
determine the effects of aberrations on these results. 

These numbers were calculated at LC = 2. If we make LC smaller 
the effect of the random lens displacements becomes less but the effect 
of correlated bends becomes larger. 7 At very small LC, the effect of 
spot size growth' becomes important. If we increase LC the effect of 
correlated bends is reduced 7 but the effect of random lens displacements 
is increased. Clearly there is an optimum LC depending on line con- 
struction tolerances and line laying tolerances. It appears this optimum 
is near LC = 2 in a typical case. 

Fig. 2 shows rms expected beam deviation as a function of n for a 
confocal line. This clearly shows how the distance between redirectors 
must be reduced as the errors in / and L become larger. 



VII. ACKNOWLEDGMENTS 



This problem and the use of Kronecker products to solve it was 
suggested by H. E. Rowe. His suggestions and comments were very 
helpful. Helpful suggestions by D. Marcuse are also acknowledged. 













1 






















































J/\ 




-& y 


















0«— - 




1 


jne cor 
L=i 


DITIONS 
METER 
















f = 1/2 METER 

Y(y 2 ^ = 2xio -2 mm 














AVERAG 


E BEND 
a = 2 


RADIUS - 

(a L 2 +ff 


5000 m 
c 2 ) 



1000 2000 3000 4000 5000 6000 7000 8000 9000 10,000 

n OR TOTAL LINE LENGTH IN METERS 



Fig. 2 — Expected value of the beam deflection as a function of transmission 
line length. 



408 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1966 

APPENDIX 

Analysis of Matrices for Bent Waveguide Axis 

We are interested in the last two terms of (5) since they contain the 
correlation between lateral lens displacements. These terms are 

((MJtn-l) X V n ) + (V » X (IA-.)). 

For simplicity let us consider only ((M n R n -i) X V n ), since the two 
terms are very similar. By repeated substitution of 



into 



R n = M n Rn-l + V n 
(MnRn-r) X Vn , 



we find 
(M n R n -i) XV n = {M n PM X V n + £ (M n P k V k -i) X V n (19) 

fc=2 

where 

P k = M n -iMn-lMn-Z ' • ' Mfc+jM* . 

If we take the expected value of (19), the first term on the right side 
is zero since it is the product of independent terms and (V n ) = 0. 

We must look at the individual elements of (M n P k V k -i) X V n . We 
can write the matrix Pk as 



P k = 



P*(l,l) P*(1.2) 



_P*(2,1) Pfc(2.2)_ 



The pk contain L„_i , L n -2 , • • • , Lk, and C„_i , C„_2 , • • • , Ck but only 
to the first power. 



M n Pk = 
and 



P*(l,l) + LnfkW.l) 

— C n Pk(l,l) + Pfc(2,l) 

(1 - L n C n ) 



P*(l,2) + I/nPfc(2,2) 

— CnPftd.S) + Pfc(2,2)(l — -W„) 



(M,P t 7w) X 7„ = 





Pft(l,2)C n + Pk(i,2)L n C n 
_~Pk(l,i)C n + P*(2,2)(w + i/nCn )_ 



Ck-lSk-lSn . 
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From this result we can write the last two terms in (5) as 

((M n R n ^) X V n ) + (V n X (MnRn-l)) 





= £ 



(pm,2))C + LC(p k tf,2)) 

(Pfc(l,2))C + LC(p k (2.2)) 
2C 2 (1 + (Tc){L(p k {2,2)) - (PH1,2))) 



C(Sk-l8 n ). (20) 



+ 2(p fc(2 ,2)>C. 

Since the pk contain the L's and C's only to the first power in each, 
(p k ) will contain only L and C and will not contain <r L 2 and a c 2 . If we 
neglect o- c as compared to 1 (the same approximation is used in (V n X 
V n )), then <(M„i? n _i) X 7»> + (7„ X (I»fl»-i)) does not contain 

2 2 

0-L Or (T C • 

We will now find under what conditions ((M n R n -i) X V n ) + (V n X 
(M„/2„_i)) is independent of n. 
The n dependence of (20) is in the terms 



S (P*(1.2))<S*-1 s n ) 



k-2 



and 



Since 



2Z (P*(2,2))(S*-1 S n ). 



<S> = <M)"- fc 

L sin (n — fc) 9 



sin 5 

L 2 C sin (n - k) d 
2 sin 



+ L cos (n — fc) 6. 



<P*(1,2)> = 
<Pfc(2.2)) = - 

And we can write 

—f— S (PHi,2))(sk-i s„) = (sis n ) sin (n — 2) 

L fc=2 

+ (s 2 s„) sin (n — 3) ■ ■ • + (s„_2S„) sin + <s„_ia„) 0. 
We will assume, as did Marcuse, that 

(s k s n ) = f(n — k). 
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The correlation depends on the distance between the lenses. We also 
assume that 

f( n - k) = for n - k > N. 

That is, the correlation length is finite and extends only N lenses away. 
This means the waveguide axis is a series of random bends, the "average 
length" of each bend is NL. Therefore, if n > N we can write 

^ S <P*(i.2))(s*-i s„> = /(2) sin 9 + /(3) sin 29 

L fc=2 

■■■ f (N - 1) sin (N - 2)9 + f{N) sin (N - 1)9, 
which is not a function of n. We can write a similar equation for 

n 

/J (Pfc(2,2))(Sfc-l S n ). 
fc=2 

However, we must consider all n down to 1. For these small 
n, ((MnRn-i) X V n ) + (V n X (MnRn-i)) will be a function of ra. This 
means that all bends contribute the same to the output beam displace- 
ment except the initial bend which is within NL of the beginning of the 
transmission line. However, if we assume that n^> N (the average bend 
length is much smaller than the length of the transmission line) the 
contribution of this initial bend will be very small and can be neglected. 

In summary, the conditions imposed on the transmission line axis 
are that it is composed of a series of random bends whose average length 
is NL. The average bend length is much smaller than the total length 
of the line. We have neglected the effect of any bend within NL of the 
beginning of the line. These are essentially the same conditions used by 
Marcuse, 7 Berreman, 6 and Unger 8 in their analyses of random bends of 
the transmission line. 

Under these conditions ((MX-i) X V„> and (V n X (M„R„-i)) 
are independent of n. 
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